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g : ABSTRACT 
On 

We examine the sensitivity of the spatial morphologies of galaxy clusters to 0,q and 
P{k) using high-resolution N-body simulations with large dynamic range. Variants 
of the standard CDM model are considered having different spatial curvatures, 
SCDM (n = 1), OCDM (fi = 0.35), LCDM (fi = 0.35, A = 0.65), and different 
normalizations, erg. We also explore critical density models with different spectral 
indices, n, of the scale-free power spectrum, P{k) oc k n . Cluster X-ray morphologies are 
quantified with power ratios (PRs), where we take for the X-ray emissivity j gas oc p 2 DM , 
which we argue is a suitable approximation for analysis of PRs. We find that Oo 
primarily influences the means of the PR distributions whereas the power spectrum (n 
and as) primarily affects their variances: log 10 (-P3/-Po) is the cleanest probe of Qq since 
its mean is very sensitive to Qq but very insensitive to P(k). The PR means easily 
distinguish the SCDM and OCDM models, while the SCDM and LCDM means show 
a more modest, but significant, difference (~ 3a). The OCDM and LCDM models 
are largely indistinguishable in terms of the PRs. Finally, we compared these models 
to a sample of ROSAT clusters and find that the PR means of the SCDM clusters 
exceed the ROSAT means with a high formal level of significance (~ 4a). Though 
the formal significance level of this p\)M I X-ray comparison should be considered 
only approximate, we argue that taking into account the hydrodynamics and cooling 
will not reconcile a discrepancy this large. The PR means of the OCDM clusters are 
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consistent, and the means of the LCDM clusters are marginally consistent, with those 
of the ROSAT clusters. Thus, we conclude that cluster morphologies strongly disfavor 
fio = 1> CDM while favoring low density, CDM models with or without a cosmological 
constant. 

Subject headings: galaxies: clusters: general - galaxies: evolution - galaxies: structure 
- X-rays: galaxies - cosmology: theory 

1. Introduction 

The quest for VLq, the current ratio of the mean mass density of the universe to the critical 
density required for closure, has been a focus of the research efforts of many astrophysicists 
involving a variety of different techniques. At present, most observational evidence suggests a 
universe with sub-critical matter density, perhaps with a cosmological constant making up the 
difference required for a critical universe (e.g., Coles & Ellis 1994; Ostriker & Steinhardt 1995). 
The possibility of measuring Q using the amount of "substructure" in galaxy clusters has thus 
generated some interest, "This is a critical area for further research, as it directly tests for Q in 
dense lumps, so both observational and theoretical studies on a careful quantitative level would be 
well rewarded." (Ostriker 1993). 

Early analytical work (e.g., Richstone, Loeb, & Turner 1992) and simulations (Evrard et al. 
1993; Mohr et al. 1995) found that the morphologies of X-ray clusters strongly favored f^o ~ 1 
over low-density universes. Along with POTENT analysis of cosmic velocity fields (e.g., Dekel 

1994) , these substructure analyses were the only indicators in support of a critical value of Slo- 
However, the analytical results (e.g., Kauffmann & White 1993; Nakamura, Hattori, & Mineshige 

1995) , simulations (e.g., Jing et al. 1994), and morphological statistics (e.g., Buote Sz Tsai 1995b) 
have been criticized rendering the previous conclusions about S7o uncertain. 

Buote & Tsai (1995b, hereafter BTa) introduced the power ratios (PRs) for quantifying the 
spatial morphologies of clusters in terms of their dynamical states. The PRs essentially measure 
the square of the ratio of a higher order moment of the two-dimensional gravitational potential 
to the monopole term computed within a circular aperture, where the radius is specified by a 
metric scale (e.g., 1 Mpc). Buote k Tsai (1996, hereafter BTb) computed PRs of ROSAT X-ray 
images for a sample of 59 clusters and discovered that the clusters are strongly correlated in PR 
space, obeying an "evolutionary track" which describes the dynamical evolution of the clusters 
(in projection). Tsai & Buote (1996, hereafter TB) studied the PRs of a small sample of clusters 
formed in the hydro dynamical simulation of Navarro, Frenk, & White (1995) and verified the 
interpretation of the "evolutionary track". In contrast to the previous studies (e.g., Richstone et 
al. 1992; Mohr et al. 1995), TB concluded that their small cluster sample, formed in a standard 
Clo = 1, CDM simulation, possessed too much substructure (as quantified by the PRs) with respect 
to the ROSAT clusters, and thus favored a lower value of SIq- 
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However, a statistically large sample of clusters is important for studies of cluster morphologies. 
The PRs are most effective at categorizing clusters into different broad morphological types; 
i.e. the distinction between equal-sized bimodals and single-component clusters is more easily 
quantified than are small deviations in ellipticities and core radii between single-component clusters 
(see BTa). The efficiency of the PRs at classifying clusters into a broad range of morphological 
types is illustrated by their success at quantitatively discriminating the ROSAT clusters along 
the lines of the morphological classes of Jones & Forman (1992) (see BTb). There is a lower 
frequency of nearly equal-sized bimodals in the ROSAT sample than clusters with more regular 
morphologies. Hence, to make most effective use of the PRs the models need to be adequately 
sampled (i.e. simulations have enough clusters) to ensure that relatively rare regions of PR-space 
are sufficiently populated. 

In this paper we build on the previous studies and investigate the ability of the PRs to 
distinguish between models having different values of S7o • Unlike the previous theoretical studies 
of cluster morphologies mentioned above, we also consider models having different power spectra, 
P(k), since P{k) should affect the structures of clusters as well. At the time we began this project 
it was too computationally costly to use hydrodynamical simulations to generate for several 
cosmological models a large, statistically robust, number of clusters with sufficient resolution. To 
satisfy the above criteria and computational feasibility we instead used pure N-body simulations. 

The organization of the paper is as follows. We discuss the selection of cosmological models 
in §2.1; the specifications of the N-body simulations in §2.2; the validity of using dark-matter-only 



simulations to generate X-ray images and the construction of the images in §2.3; and computation 



of the PRs in §2.4. We analyze the models having different values of Slo and a cosmological 
constant in jjj^, and models with different spectral slopes and as in §||. The implications of the 
results for all of the models and comparison of the simulations to the ROSAT sample of BTb is 
discussed in §||. Finally, in §|6]we present our conclusions. 



2. Simulations 

2.1. Cosmological Models 

To test the sensitivity of cluster morphologies to the cosmological density parameter due to 
matter, Qq, and the power spectrum of density fluctuations, P(k), we examined several variants of 
the standard Cold Dark Matter (CDM) model (e.g., Ostriker 1993). In Table | we list the models 
and their relevant parameters: Oo; Ao = A/3-Hq, where A is a cosmological constant and Hq is the 
present value of the Hubble parameter; the spectral index, n, of the scale-free power spectrum of 
density fluctuations, P (k) oc k n ; and erg, the present rms density fluctuations in spheres of radius 
8/i -1 Mpc, where h is defined by Hq = lOO/i km s _1 Mpc -1 . 

The parameters of the open CDM model (OCDM) and low-density, flat model (LCDM) 
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were chosen to be consistent with current observations (e.g., Ostriker & Steinhardt 1995). Their 
normalizations were set according to the erg — flo relationship of Eke, Cole, & Frenk (1996) to 
agree with the observed abundance of X-ray clusters. The biased CDM model (BCDM) was also 
normalized in this way. However, the BCDM simulation, because it has £Iq = 1, necessarily has 
poorer resolution (i.e. fewer particles per cluster) than the OCDM and LCDM models due to 
the fixed box size of our simulations (see §[2.3]) . For the purposes of our investigation of cluster 
morphologies it is paramount to compare simulations having similar resolution. Hence we use the 
SCDM model (with a$ = 1) as our primary f^o = 1 simulation for analysis, which has resolution 
equivalent to the OCDM and LCDM simulations. (We show in § [4.4| that the means of the PR 
distributions for BCDM and SCDM are very similar which turns out to be most important for 
examining the effects of rZo-) Hence, the SCDM, OCDM, and LCDM models allow us to explore 
the effects of f2o and Ao on the cluster morphologies; comparing SCDM and BCDM provides 
information on the influence of <7g. 

We explore the effects of different P{k) on the PRs using the scale-free models, which have 
different n from SCDM. For the scale-free models we normalized each to the same characteristic 
mass, M*, defined to be (Cole & Lacey 1996) the mass scale when the linear rms density 
fluctuation is equal to <5 C , the critical density for a uniform spherically symmetric perturbation to 
collapse to a singularity. For Qq = 1, the linear theory predicts 5 C ~ 1.686 (e.g., Padmanabhan 
1993). We take the SCDM model with as = 1 as a reference for these scale- free models which 
gives a characteristic mass of 10 14 Mq. This procedure allows a consistent means to normalize the 
scale-free models relative to each other on the mass scales of clusters. Unfortunately, as a result 
of this normalization procedure, at earlier times the models have different large-scale power and 
thus the cluster mass functions are different for each of the models. The scale-free model with 
n = —1.5 is similar to the SCDM model and will be used to "calibrate" the scale-free models with 
respect to the other models (see Table [j]). 



2.2. N-body Cluster Sample 

We use the Tree-Particle-Mesh (TPM) N-body code (Xu 1995b) to simulate the dissipationless 
formation of structure in a universe filled with cold dark matter. The simulations consist of 128 3 
particles in a square box of width 200/i~ 1 Mpc. The gravitational softening length is 25/i~ 1 kpc 
which translates to a nominal resolution of ~ 50/i _1 kpc. This resolution is sufficient for exploring 
the structure of clusters with PRs in apertures of radii R ap ;> 0.5 Mpc; for a discussion of the 
related effects of resolution on the performance of PRs on ROSAT X-ray images see Buote &: Tsai 
(1995b, §4). All of the realizations have the same initial random phase. 

For each simulation we located the 39 most massive clusters using a version of the DENMAX 
algorithm (Bertschinger & Gelb 1991) modified by Xu (1995a). This convenient selection criterion 
yields well defined samples for each simulation and allows consistent statistical comparison 
between different simulations which is the principal goal of our present investigation. For the 
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various cosmological models we explore (see Table [l]) these clusters generally have masses ranging 
from (0.3 — 3) x h 10 Mq, which correspond to typical cluster masses observed in X-ray (e.g., 
Edge et al. 1990; David et. al. 1993) and optically (e.g., Carlberg et. al. 1995) selected samples. 

2.3. X-ray Images 

2.3.1. Motivation for j g oc p 2 D ^[ 

By letting the gas density trace the dark matter density (p ga s oc pdm) and by assuming that 
the plasma emissivity of the gas is constant, we computed the X-ray emissivity of the clusters, 
j g oc p 2 DM - Given its importance on the results presented in this paper, here we discuss at some 
length the suitability of this approximation. (Cooling of the gas is discussed in § |5.1| .) 

For clusters in the process of formation or merging, the gas can have hot spots appearing 
where gas is being shock heated (e.g., Frenk, Evrard, Summers, & White 1996). One effect of such 
temperature fluctuations on the intrinsic X-ray emissivity is that the intrinsic plasma emissivity 
will vary substantially over the cluster thus rendering j g oc p\,w a poor approximation. However, 
the intrinsic X-ray emissivity is not observed, but rather that which is convolved with the spectral 
response of the detector. For ROSAT observations of clusters with the PSPC, the plasma 
emissivity is nearly constant over the relevant ranges of temperatures (NRA 91-OSSA-3, Appendix 
F, ROSAT Mission Description), and thus temperature fluctuations contribute negligibly to 
variations in the emissivity; for previous discussions of this issue for the PRs see BTa and TB. 

A more serious issue is whether the shocking gas invalidates the p gas oc pdm approximation, 
in which case the dynamical state inferred from the gas would not reflect that of the underlying 
mass. TB, who analyzed the hydro dynamical simulation of Navarro et al. (1995a), showed that 
the PRs computed for both the gas and the dark matter gave similar indications of the dynamical 
states of the simulated clusters (see §4 of TB). In particular, this applied at early times when the 
clusters underwent mergers with massive subclusters. Hence, our approximation for the X-ray 
emissivity should be reasonable even during the early, formative stages of clusters. 

Another possible concern with setting p gas oc pdm is that a gas in hydrostatic equilibrium, 
which should be a more appropriate description for clusters in the later stages of their evolution, 
traces the shape of the potential of the gravitating matter which is necessarily rounder than the 
underlying mass; if the gas is rounder, then the PRs will be smaller. However, the core radius (or 
scale length) of the radial profile of the gas also influences the PRs. In fact, clusters with larger 
core radii have larger PRs; see BTa who computed PRs for toy X-ray cluster models having a 
variety of ellipticities and core radii. 



2 See Buote & Tsai (1995a) for a related discussion of the evolution of the shape of the gas and dark matter in the 
Katz & White (1993) simulation. 
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When isothermal gas, which is a good approximation for a nearly relaxed cluster, is added 
to the potential generated by an average cluster formed in a flo = 1, CDM simulation, the gas 
necessarily has a larger core radius than that of the dark matter (see Figure 14 of Navarro, Frenk, 
& White 1995b). Hence, a gas in hydrostatic equilibrium will have a larger core radius than that 
of the dark matter, at least in the context of the CDM models we are studying. Considering the 
competing effects of smaller ellipticity and larger core radii (factors of 2-3 in each), and from 
consulting Table 6 of BTa, we conclude that no clear bias in the PRs is to be expected by assuming 
that the gas follows the dark matter. In further support of this conclusion are the similarities 
of the morphologies of clusters in the N-body study of Jing et al. (1995). Using centroid-shifts 
and axial ratios they find similar results when p gas oc pdm and when the gas is in hydrostatic 
equilibrium (see their Figures 5, 6, and 8). 

2.3.2. Construction of the Images 

Having chosen our representation for the X-ray emissivity, we then generated two-dimensional 
"images" for each cluster. A rectangular box of dimensions 4 x 4 x 10 /i -3 Mpc 3 with random 
orientation was constructed about each cluster. We converted the particle distribution for each 
cluster to a mass density field using the interpolation technique employed in Smoothed Particle 
Hydrodynamics (SPH) (e.g., Hernquist & Katz 1989), from which the X-ray emissivity was 
generated, j g cx p 2 DM . The SPH interpolation calculates the density at a grid point by searching 
for the nearest neighbors and is thus more robust and physical than other linear interpolation 
schemes like Cloud-in-Cell. For our SPH interpolation we use 20 neighbors and the spline kernel 
described in Hernquist & Katz. The interpolation result is independent of the cell we choose for 
the X-ray emissivity calculations. 

Typically, the boxes contained ~ 2500 particles for a cluster; e.g., SCDM (1799-3965), OCDM 
(1359-3656), LCDM (1439-3877), and BCDM (603-1740). We projected the emissivity along the 
long edge of the box into a square 4x4 h~ 2 Mpc 2 "image" consisting of "pixels" of width 20fo -1 
kpc. This pixel width was chosen sufficiently small so as not to inhibit reliable computation the 
PRs. 

We do not add statistical noise or other effects associated with real observations to the X-ray 
images since our principal objective is to examine the intrinsic response of cluster morphologies to 
different cosmological parameters. However, the investigation of observational effects on the PRs 
by BTa, and the derived error bars on the PRs from ROSAT clusters by BTb, do not show any 
large systematic biases; comparison of the simulations to the ROSAT cluster sample is discussed 
in §1 
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2.4. Power Ratios 

The PRs are derived from the multipole expansion of the two-dimensional gravitational 
potential, ty(R,(f)), generated by the mass density, E(R,(f>), interior to R, 

( 1 \ 00 1 

^>(R,(f)) = -2Ga ln — ) - 2G 2J — ^7 (a m cos + b m sin m</>) , (1) 



R J , mi? 

771=1 

where is the azimuthal angle, G is the gravitational constant and, 

a m {R) = f E(x') (R') m cos m<j)'d 2 x', (2) 
b m (R) = [ ^{x')(R') m smm<t)'d 2 x'. (3) 

JR'<R 

Because of various advantageous properties of X-ray images of clusters, we associate the surface 
mass density, E, with X-ray surface brightness, Ex (which is derived from the projection of pjjM 
- see previous section); for more complete discussions of this association see BTa and TB. The 
square of each term on the right hand side of eq. (|l|) integrated over the boundary of a circular 
aperture of radius R ap is given by (ignoring factors of 2G), 

Pm = 2m 2 R 2 ™ + 6 ™) ' (4) 

for m > and, 

P = [a ln{R ap )} 2 , (5) 

for m = 0. It is more useful for studies of cluster structure to consider the ratios of the higher 
order terms to the monopole term, P m /Po, which we call "power ratios" (PRs). By dividing each 
term by Pq we normalize to the flux within R ap - Because for clusters, P m /Po <ti 1 for m > (e.g., 
BTb), it is preferable to take the logarithm of the PRs, 

Pm 

PR m = log 10 -^, (6) 
"0 

which we shall henceforward analyze in this paper. 

Since the P m depend on the origin of the chosen coordinate system, we consider two choices 
for the origin. First, we take the aperture to lie at the centroid of Sx; i.e. where Pi vanishes. 
Of these centroided PR m , PR2, PR3, and PR4 prove to be the most useful for studying cluster 
morphologies (see BTa). In order to extract information from the dipole term, we also consider the 
origin located at the peak of Ex- We denote this dipole ratio by P^^ / pjf k \ and its logarithm 

(vk) 

PR\ , to distinguish it from the centroided power ratios. 

To obtain the centroid of Ex in a consistent manner for all clusters we adopted the following 
procedure. First, when projecting the cluster (see § |2.3| ), the cluster was roughly centered on 
the X-ray image by eye. For each image we computed the centroid in a circular aperture with 



- 8 - 



R ap = 1.5 h Mpc located about the field center. This centroid was then used as our initial 
center for each cluster; see BTa for a description of how the peak of £x is located. 

In addition to considering the PR m individually, we also analyze the cluster distributions 
along the "evolutionary tracks" in the (PR2, PRa) and (PR2, PR3) planes obeyed by the ROSAT 
clusters of BTb. We refer the reader to TB for a detailed discussion of the cluster properties along 
the evolutionary tracks. 

Using the augmented Edge et al. (1990) sample of BTb we recomputed the lines defining 
the evolutionary tracks of the ROSAT data for the l/igQ 1 Mpc apertures. This was done since 
TB selected a subset of the clusters based on PR m measurement uncertainty rather than flux. 
Following TB we fit PR4 = a + bPR2 considering the uncertainties in both similar fit 

was done for the (PR2, PR3) plane. We obtained a = —0.92, b = 1.18 for the (PR2, PRa) track 
which we denote by PR2-A- Similarly, for the (PR2, PR3) track, which we denote by PR2-3, we 
obtained a = —0.49, b = 1.16. These results are nearly the same found by TB for the slightly 
different sample. 

To facilitate comparison to a previous study of the PR m of ROSAT clusters (BTb), we 
compute PR m of the simulated clusters in apertures ranging in radius from 0.5/i^Q 1 Mpc to 1.5/igQ 1 
Mpc (H = 80/i 80 km s" 1 Mpc" 1 ) in steps of 0.25h^ Mpc; i.e. (0.4 - l.2)h~ l Mpc in steps of 
0.2/j" 1 Mpc. We refer the reader to §2 of BTa and §2 of TB for discussions of the advantages of 
using a series of fixed metric aperture sizes to study cluster morphologies. 



3. PR m for Models with Different J)q an d Ao 

First we consider clusters formed in the SCDM, OCDM, and LCDM models. Contour plots 
for 16 of the clusters formed in each of the models are displayed in Figures |], and |3[ In 
Figure ^| we show the (PR2, PRa) plane for the 0.5/i^q 1 Mpc and l-0hg O Mpc apertures; the 
SCDM model appears in each plot for comparison. The clusters in each of the models exhibit 
tight correlations very similar to the evolutionary tracks of the ROSAT clusters (BTb) and 
the simulated hydrodynamic clusters (f2o = 1) of Navarro et al. (1995a) studied by TB. Along 
the evolutionary tracks a shift in the means of the PRm is easily noticeable in the 0.5/igQ 1 Mpc 
aperture, being most apparent for the SCDM-OCDM models. The spread of the PR m along the 
track in the lM^ Mpc aperture for SCDM-OCDM also appears to be different. The distributions 
perpendicular to PR2-A do not show discrepancies obvious to the eye. 

At this time we shift our focus away from the evolutionary tracks and instead analyze 
the individual PR m distributions, which prove to be more powerful for distinguishing between 
cosmological models as we show below. We give the individual PR m distributions of clusters in the 
three models for the 0.5h$Q Mpc and 1.0hg Mpc apertures in Figures || and || We found it most 
useful to compare these distributions in terms of their means, variances, and Kolmogorov-Smirnov 
(KS) statistics. For the number of clusters in each of our simulations (39) higher order statistics 
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like the skewness and kurtosis are unreliable "high variance" distribution shape estimates (e.g., 
Bird & Beers 1993). We did consider more robust statistics like the "Asymmetry Index" (AI), 
which measures a quantity similar to the skewness, and the "Tail Index" (TI), which is similar 
to the Kurtosis (Bird & Beers 1993). However, we found that they did not clearly provide useful 
information in addition to the lower order statistics and KS test, and thus we do not discuss them 
further. ^\ 

The means and standard deviations for the (0.5,0.75, 1.0) /i^j Mpc apertures are listed in 
Table As a possible aid to understanding the relationships between the values in Table 0, 
we plot in Figure [?] the standard deviation vs the mean for the 0.75/igQ 1 Mpc aperture. We do 
not present the results for the larger apertures because they did not significantly improve the 
ability to distinguish between the models. Moreover, we found that PR^ and PR^ k ^ do not 
provide much useful information in addition to PR2 and PR3. Generally P-R4 tracks the behavior 
of PR2, though showing less power to discriminate between models; the similarity to PR2 is 
understandable given the strong correlation shown in Figure ||. Likewise, PR^f^ is similar to, but 
not quite so effective as, PR3. For compactness, thus, we shall henceforward mostly restrict our 
discussion to results for PR2 and PR3 in the (0.5,0.75, 1.0)h$ Q Mpc apertures. 

We compare the means, standard deviations, and total distributions of the models in Table 
H using standard non-parametric tests as described in Press et al. (1994). The Student's t-test, 
which compares the means of two distributions, computes a value, pt, indicating the probability 
that the distributions have significantly different means. Similarly, the F-test, which compares 
the variances of two distributions, computes a value, pp, indicating the probability that the 
distributions have significantly different variances. Finally, the KS test, which compares the 
overall shape of two distributions, computes a value, PkSi indicating the probability that the 
distributions originate from the same parent population; the probabilities listed in Table [3| are 
given as percents; i.e. decimal probability times 100. Note that for the cases where the F-test 
gives a probability less than 5% we use the variant of the t-test appropriate for distributions with 
significantly different variances (i.e. program tutest in Press et. al.). 

3.1. SCDM vs. OCDM 

As is clear from inspection of Figures I - 0, and Tables g - §, the means of the PR m of the 
SCDM model exceed those of OCDM. In terms of the t-test the significance of the differences is 
very high. Of all the PR m , generally the means of PR2 and PR3 exhibit the largest significant 
differences; the most significant differences are seen for PR3 in the 0.5/i^q 1 Mpc aperture, 
Pt = 0.02%, and for PR2 in the 0.75/igQ 1 Mpc aperture, pt = 0.06%. Hence, though different in 



3 Actually, the KS test turns out not to provide much additional information over the t-test and F-test, but we 
include it for ease of comparison to previous studies; e.g., Jing et al. (1995); Mohr et al. (1995); TB. 
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all the apertures, the discrepancy in the means is most significant for the smallest apertures, 
(0.5, 0.75)k^ Mpc. 

The variances of PR3 in the SCDM model are essentially consistent at all radii with their 
corresponding values in OCDM. However, for PR2 the variances are consistent at 0.75/i^q 1 Mpc, 
but marginally inconsistent at (0.5, 1.0) Mpc (and inconsistent at (1.25, 1.5)h^ Mpc). 

The KS test generally indicates a significant difference in SCDM and OCDM when also 
indicated by the t-test, or the t-test and F-test together. The level of discrepancy is usually not 
as significant as given by the t-test, except when pf is small as well. Since the KS test does not 
indicate discrepancy when both the t-test and F-test indicate similarity, we conclude that higher 
order properties of the PR distributions are probably not very important for the SCDM and 
OCDM models (at least for the samples of 39 clusters in our simulations). Since this qualitative 
behavior holds for the other model comparisons, we shall not emphasize the KS tests henceforward 

Finally, in terms of the various significance tests we find that the PR2-4 distribution 
essentially gives a weighted probability of the individual PR2 and PR4 distributions; i.e. it 
does not enhance the discrepancy in the individual distributions. Perpendicular to PR2-A the 
distributions are consistent. The same behavior is seen for PR2-3 as well. This behavior is seen 
for the remaining model comparisons in this section so we will not discuss the joint distributions 
further. 



3.2. SCDM vs. LCDM 

The PR m means for the SCDM clusters also systematically exceed those in the LCDM model, 
however the significance of the difference is not as large as with the OCDM clusters. The largest 
discrepancy is observed for PR2 in the (0.75, 1.0)/i 80 1 Mpc apertures for which pt = (0.7%, 0.8%). 
The other PR m show only a marginal discrepancy in the means. For apertures (0.5, 0.75)hgQ Mpc, 
PR3 has pt = (4%, 3%), but is quite consistent at larger radii. The variances for the SCDM and 
LCDM models are consistent for essentially all radii and all PR m . 

3.3. OCDM vs. LCDM 

The PR m means for the LCDM clusters appear to systematically exceed those in the OCDM 
model, however the formal significances of the differences are quite low. The means are entirely 
consistent at all radii for PR2. However, PR3 shows a marginal difference in the 0.5/igQ 1 Mpc 
aperture (jpt = 9%). The variances of the PR m of the OCDM and LCDM models behave similarly 
as with the SCDM and OCDM comparison above, as expected since the SCDM-LCDM variances 
are essentially identical. However, the degree of discrepancy is not as pronounced. 
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3.4. Performance Evaluation I. 

The means of the individual PR m distributions generally exhibit the most significant 
differences between the SCDM, OCDM, and LCDM models; the variances are much less sensitive 
to the models, with PR3 showing no significant variance differences. The larger means for the 
PR m in the SCDM models are expected from the arguments of, e.g., Richstone et al. (1992). That 
is, in a sub-critical universe the growth of density fluctuations ceased at an early epoch and so 
present-day clusters should show less "substructure" than in an f^o = 1 universe where formation 
continues to the present. Clusters with more structure will have systematically larger values of the 

PRm- 

The PR m whose means show the most significant differences between the models are PR2 
and PR3, where PR2 typically performs best for apertures (0.75, 1.0)/tg Mpc and PR3 is most 
effective for (0.5, 0.7b)hjQ Mpc. Although useful, PpJf k ^ is often the least effective PR m for 
differentiating models in terms of its mean; this relatively weak performance of the dipole ratio 
with respect to other moments is echoed in the results of Jing et al. (1995) who found that their 
measure of an axial ratio performed better than a centroid shift for discriminating between models 
(see their tables 3-6). 



4. PR m for Models with Different n and eg 

In this section we investigate the effects of different power spectra for models otherwise 
conforming to the specifications of the SCDM model. First, we examine models with different 
spectral indices of the scale-free power spectrum (P(k) oc k n ), n = 0, —1, —1.5, —2. Then we 
examine the BCDM model which has a lower power-spectrum normalization as expressed by as- 
As in the previous section, we find the (0.5, 0.75, 1.0) h^Q Mpc apertures to be more useful than the 

(pk) 

larger apertures, and that PR<± and PR\ do not provide much useful information in addition to 
that provided by PR2 and PR3. Hence, for compactness we again mostly restrict the discussion 
to PR2 and PR3 in the smaller apertures. 



4.1. n = -15 vs. SCDM 

Before analyzing the PR m of models with different n we calibrate the scale-free models by 
comparing the n = —1.5 scale-free model to the SCDM model since they should have similar 
properties (see §|2.1|). We find that the means, variances, and KS statistics of the centroided PR m 
for the SCDM and n = —1.5 models are entirely consistent for all aperture radii with only one 
possible exception. The variances of PR3 exhibit a marginal (jpp = 5%) discrepancy in the 0.5/igg 1 
Mpc aperture. The significance of this variance discrepancy should be treated with caution given 
the complete consistency of the means (pt = 31%) and KS (32%) test at this radius as well as the 
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consistency of all the tests at all the other radii investigated. Hence, the cluster morphologies of 
the SCDM and n = —1.5 models are very consistent expressed in terms of the centroided PR m 
(m = 2,3,4). 



4.2. n = vs. n = — 2 

In Figure ^| we plot for the n = 0, —2 models the PR correlations for m = (2, 3) and m = (2, 4) 
in the (0.5, l.O)/^ 1 Mpc apertures. Histograms for the individual PR m in these apertures are 
displayed in Figures | and 0. Table | lists the means, Figure [ll| plots the standard deviations 
versus the averages of the PR m in the 0.75/igQ 1 Mpc aperture, and Table |3| gives the results of the 
significance tests. 

The means of PR3 are very consistent for the n = 0, —2 models at all radii examined. Those 
for Pi?2 may show some differences in their means, with the n = —2 models perhaps having 
systematically smaller values. The significance of the different means for PR2 is only formally 
marginal, with p t = (11%, 4%, 10%)% for aperture radii (0.5,0.75, 1.0)K^ Mpc. 

However, the possible small differences in the means of PR2 are dwarfed by the corresponding 
highly significant differences in its variances. Generally the variances for all the PR m in all 
the apertures are smaller for the n = —2 clusters. The most significant variance differences are 
observed for PR2 which has pf < 1% in (0.75, l.O)/^ 1 Mpc apertures and pf ~ 5% for 0.5/igg 1 
Mpc . The variances for PR3 show differences but at a lower level of significance and only in the 
(0.5,0.75, l.Q)K$ Mpc apertures; i.e. p t = (3%,1%,4%). 

Similar to what we found in §|3[ the differences implied by the KS test generally follow the 
significances implied by the t-test and F-test; i.e. higher order effects in the distributions are 
probably not overly important (at least for our sample sizes of 39 clusters). Moreover, again we 
find that analysis of the PR m in terms of the evolutionary tracks does not add useful information 
to the previous results. The mean and variance effects for the individual PR m translate to very 
similar behavior along Pi?2-4- The direction perpendicular to PR2-4 is essentially consistent for 
all of the tests. As a result, we do not emphasize the KS tests or the evolutionary tracks further. 



4.3. Intermediate n 

The behavior for other n is similar, but depends to some extent on the range examined. We 
find that the range of n which accentuates differences in the PR m is between n = 0,-1. Over the 
range n = 0,-1 the discrepancy of means for PR2 essentially follows that of the full n = 0,-2 



discussed in §4.2. However, the variances are not so highly discrepant as before, with pp = 3% for 
PR2 for aperture radii (0.5, 0.75)/ig 1 Mpc; elsewhere the variances of PR2 are consistent between 
the n = 0, — 1 models. Over the n = 0, — 1 range PR3 is consistent for all statistics at all radii 
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examined. 

The PR m exhibit very few differences over the range of indices n = —1,-2. For all radii 
examined the means and KS statistics are consistent for all the PR m - However, the variances do 
show some marginal differences. The 1.0/igQ 1 Mpc apertures has the most significance difference 
where pp = 1.5% for both PR2 and PR3. Also in the 0.75/jgQ 1 Mpc aperture PR2 has pp = 7%. 
Otherwise the variances of these PR m ar e consistent. (We mention that the variance of PR2 for 
the n = —1.5 model in the 0.75/igQ 1 Mpc aperture lies above that for the n = — 1 model in Figure 



11, but the difference is not statistically significant.) 



4.4. SCDM vs. BCDM 

Now we consider the £Iq = 1, CDM model with a lower normalization, as = 0.51, which we 
refer to as the biased CDM model, BCDM. The means and variances for the BCDM model are 
listed in Table [2], the standard deviation versus the average PR m in the 0.75/igQ 1 Mpc aperture are 
plotted in Figure [ll], and the results for the significance tests in comparison to SCDM are given in 
Table |. 

The means of all the PR m at all aperture radii are consistent for the SCDM and BCDM 
models. The PR m variances of the SCDM clusters generally exceed those of the BCDM clusters. 
The significance levels of the differences are only marginal {pp ~ 3%) and appear to be most 
important in the 0.75/igQ 1 Mpc aperture. 

It is possible that the slight variance differences between the SCDM and BCDM models are 
due to the difference in resolution between the two simulations; i.e. the clusters in the BCDM 
simulations contain about half the number of particles of the SCDM clusters. We would expect 
that the effects of resolution would be most important in the smallest apertures (which we do 
observe), although we would probably expect that the means as well as the variances would 
be affected (which we do not observe). We mention that the BCDM model performs virtually 
identically to the SCDM model when compared to the OCDM and LCDM models. 



4.5. Performance Evaluation II. 

The variances of the PR m show the most significant differences between models with different 
power spectra; PR2 generally has the most sensitive variances over the parameter ranges explored. 
Decreasing n and ag both decrease the PR m variances, the differences being of similar magnitude 
for the n = 0, —1 models and the SCDM and BCDM models. The means of the PR m are much 
less sensitive to the models with different n and erg, with PR2 showing the largest significant 
differences which are always less than differences in the variances. No significant differences in the 
means are observed for PR3 over the range of power spectra studied. 
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The predominant effect of the power spectrum on the variances of the PR m is intriguing. It is 
reasonable that when the amount of small-scale structures is reduced (smaller n) or the population 
of cluster-sized structures is made more uniform (smaller a%) that the PR m distributions would 
also be more uniform. The observed low sensitivity of the PR m means to the power spectra is also 
reasonable since on average the PR m means should only be affected by the rate of mass accretion 
through the aperture of radius R ap , not by the sizes of the individual accreting clumps. 

5. Discussion 

In the previous sections we have seen that differences in Qq and P{k) in CDM models are 
reflected in the spatial morphologies of clusters when expressed in terms of the PR m - For the 
purposes of probing Oo> our analysis indicates that PR3 is the best PR since its mean is quite 
sensitive to f2 but very insensitive to P(k). It is advantageous to also consider PR2 when a 
cosmological constant is introduced since its means differ for the SCDM and LCDM models by 
~ 3a whereas PR3 only distinguishes the models at the ~ 2a level. The marginal dependence 
of the mean of PR2 on P(k) is not overly serious for studying differences in because the 
differences in means due to P{k) are always accompanied by larger, more significant differences in 
the variances; i.e. different means for PR2 but consistent variances should reflect differences only 
in Qo- The best apertures for segregating models are generally (0.5,0.75, 1.0)^ Mpc. 

A few previous studies have examined the influence of an d Ao on the morphologies of 
galaxy clusters. Perhaps the most thorough investigation is that of Jing et. al. (1995) who used 
N-body simulations of a variety of CDM models, including versions similar to our SCDM, OCDM, 
and LCDM, to study variations of center-shifts and axial ratios. Jing et. al. reached the same 
qualitative conclusions as we do; i.e. the SCDM model is easily distinguished from OCDM and 
LCDM because it produces clusters with much more irregular morphologies than than the others. 
However, Jing et al. obtained infinitesimal KS probabilities for the axial ratio when comparing 
SCDM to OCDM and LCDM, a level of significance orders of magnitude different from that found 
in this paper. The source of this discrepancy is unclear given the qualitative similarities of their 
axial ratio and our PR2- The disagreement may arise from differences in numerical modeling 
between the simulations; i.e. the results of Jing et al. are derived from simulations with a larger 
force resolution (O.l/i -1 Mpc), and smaller particle number for the non-SCDM models (64 3 ) than 
in our simulations, and have clusters which visually do not show the rich structures seen in our 
simulations (Figures [l], [2], and ||) . 

The qualitative results of Jing et al. agree with the hydrodynamic simulations of Mohr et 
al. (1995) who also used center shifts and axial ratios as diagnostics for "substructure". If we 
visually estimate the means of the center shifts and axial ratios from Figures 6 and 8 of Jing et 
al. for their SCDM, OCDM, and LCDM models (actually OCDM with ft = 0.2 and LCDM with 
f^o = 0.2, Ao = 0.8), we find that they agree quite well with the corresponding values in Table 3 
of Mohr et al.; i.e. the results from the N-body and hydrodynamic simulations are very similar, 
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despite the many other differences between the simulations (e.g., large number of baryons in 
OCDM clusters for Mohr et al.). 

We can make a similar comparison of the PR m derived in this paper with the results from 
TB who analyzed the small sample of SCDM clusters formed in the hydrodynamic simulation of 
Navarro et al. (1995a). We find that the means (and variances) of the PR m computed in this 
paper are very similar to those of the hydrodynamic clusters; e.g., the mean for Pi?2-4 f° r l^o 1 
Mpc may be read off Figure 7 of TB which shows excellent agreement with the SCDM value we 
obtain from the N-body simulations (average PR2-A = 3.76). The quantitative similarity between 
the results, particularly between the means of the morphological statistics, for the N-body and 
hydrodynamic simulations of (Jing et al.,Mohr et al.) and (this paper, TB) suggest that it is useful 
to compare the PR m derived from N-body simulations directly to the X-ray data. 



5.1. Comparison to ROSAT Clusters 

Among the biases that need to be considered in such a comparison are the effects of cooling 
flows (e.g., Fabian 1994), selection, and noise. Cooling flows increase the X-ray emission in the 
cluster center, which has the effect on the PR m of essentially decreasing the core size of the 
cluster. Judging by the observed core radii of "regular" X-ray clusters we would expect at most 
a factor of ~ 2 difference in core radii (e.g., A401 vs. A2029 in Buote & Canizares 1996; also see 
Jones & Forman 1984). Q Changing the core radius by a factor of 2 typically changes PR2 (for 
example) by a small fraction of a decade (see Table 6 of BTa); this behavior, as we show below, 
is confirmed using a more thorough treatment. The issue of biases between X-ray-selected and 
mass-selected samples needs to be addressed with hydrodynamical simulations. The estimated 
uncertainties of the PR m for the ROSAT cluster sample of BTb, which take into account noise 
and unresolved sources, do not show any clear biases. 

In Figure ^ we display the correlations of the centroided PRm for the ROSAT sample of 
BTb in the (0.5, 1.0)hg Mpc apertures; the SCDM clusters are also plotted for a comparison. 



PR m histograms for these apertures are shown in Figure 13 and 14, along with those for the 
SCDM and OCDM models. We list the means and variances for the ROSAT clusters in Table 
||; we plot in Figure [l5| the standard deviations versus the means for the ROSAT clusters and 
models in the 0.5/i^q 1 Mpc aperture; the results of the significance tests between the ROSAT 
clusters and model clusters are given in Table [5|. We analyze the ROSAT clusters corresponding 
to the "updated Edge et. al. (1990)" flux-limited sample in BTb which gives 37 and 27 clusters 
respectively for the (0.5, l.O)/^ 1 Mpc apertures; note that all the qualitative features of the results 
we obtain below are reproduced when all of the clusters studied in BTb are used (i.e. 59 and 44 



4 Large cooling flows only appear in clusters with regular morphologies (e.g., Jones & Forman 1992; Fabian 1994; 
BTb). 
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clusters respectively). 

The means of the SCDM clusters exceed those of the ROSAT sample to a high level of 
significance, with the differences being most pronounced in the 0.5/igQ 1 Mpc aperture. The most 
significant discrepancy is for PR3 in the 0.5/i^q 1 Mpc aperture for which pt = 1.5 x 10~ 4 %. The 
variances for all the PR m except PR3 are also significantly different, with the variances of the 
SCDM clusters exceeding those of the ROSAT clusters. The SCDM model has as = 1 which is 
too high to fit other observations (e.g., Ostriker & Steinhardt 1995). The BCDM model, which has 
as = 0.51, does have PR m variances in better agreement with the ROSAT sample. However, the 
means are in essentially the same level of disagreement. In fact, PR2 has a much more significant 
mean discrepancy (pt = 1.0 x 10~ 4 %) in the O.bhgQ Mpc aperture. 

In contrast, the PR m have means that are entirely consistent for the OCDM and ROSAT 
clusters in both apertures. The variances of the centroided PR m are significantly discrepant, 
particularly in the 1.0/igQ 1 Mpc aperture, where the OCDM variances exceed the ROSAT 
variances. This suggests a lower as or n is needed to bring the variances of the OCDM models 
into agreement with the ROSAT sample. 

The PR m means of the LCDM clusters systematically exceed the ROSAT means, but 
at a lower level of significance than does SCDM. The discrepancies are only significant in the 
0.5/igQ 1 Mpc aperture, where PR3 (pt = 0.4%) and PR^ (p t = 0.2%) show the most significant 
discrepancies; the even PRm show at best a marginal discrepancy in their means (p t = 10% — 15%). 
The variances for the even PR m are also significantly different, though only in the 0.5/igQ 1 Mpc 
aperture as well. As the LCDM and SCDM variances are very similar, we expect that the variance 
differences can be largely obviated with a lower value of as- 

The difference in the means of PR 3 for the LCDM and ROSAT clusters in the 0.5/igQ 1 Mpc 
aperture, though formally significant at better than the 3a level, represents a shift of about 
one-half a decade in PR3; also, when using all 59 clusters of BTb the significance is only p t = 4% 
(~ 2a). As we have discussed earlier, it is difficult to completely account for a discrepancy of this 
magnitude by invoking, e.g., the unsuitability of the p gas oc pdm approximation, observational 
noise, or cooling flows. 

We may make a more precise estimate of the effects of cooling flows on the PRm- The 
ROSAT clusters in the augmented Edge sample all have estimated mass-flow rates (Fabian 1994) 
from which we may compute a luminosity (bolometric) due to the cooling flow following Edge 
(1989), L coo i = 3.0 x lO 41 /^ 2 MT erg/s, where M is in Mq /year and T is in keV. Comparing this 
cooling luminosity to the total cluster luminosity, L^, using the results of David et al. (1993) 
allows us to in effect remove the cooling gas from the ROSAT PR m - To a first approximation the 
cooling flow affects only Pq because the cooling emission is weighted heavily towards the aperture 
center. Hence, to approximately remove the effects of the cooling flows from the ROSAT clusters 
we reduce Po f° r each cluster by (1 — L coo i/Lb i). We find that the PR m of the ROSAT clusters 
are modified minimally, the effect being that the means of the PR m are increased by 1/10 of a 
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decade: means for PR2 and PR3 are -5.60 and -7.52 respectively in the 0.5/ig aperture; the 
variances show no significant systematic effect. These small mean shifts do reduce the significance 
of the LCDM-ROSAT discrepancy, but the discrepancy is still significant at the ~ 3<r level; e.g., 
Pt = 1.6% for PR3 and pt = 1% for PR^f k ^ in the 0.5/igQ 1 aperture, and pt = 34% for the even 

PRm- 

Although cooling flows alone cannot completely account for the differences in the ROSAT 
clusters and the LCDM model, it is very possible that when combined with the the other effects 
mentioned above a sizeable fraction of the half-decade difference could be made up which would 
in any event reduce the significance level of the difference. As a result, we believe the discrepancy 
of the LCDM- ROSAT means must be considered preliminary and await confirmation from 
appropriate hydrodynamical simulations 

On the other hand, the means of PR 3 for the SCDM and BCDM models exceed the ROSAT 
means by almost a full decade to a higher formal significance level (~ 4cr), which in light of the 
previous discussion should be considered robust. We conclude that the Oo = 1, CDM models 
cannot produce the observed PR m of the ROSAT clusters, and that the discrepancy in PR m 
means is due to being too large. This agrees with our conclusions obtained in TB for the small 
sample of clusters drawn from the hydrodynamic simulation of Navarro et al. (1995a). [] 

Our conclusions are opposite those of Mohr et al. (1995) who instead concluded that their 
Einstein cluster sample favored SCDM over both OCDM and LCDM. Given the qualitative 
agreement discussed above between the Jing et al. and Mohr et al. simulations, as well as 
between our present simulations and TB, it would seem that the discrepancy lies not in the details 
of the individual simulations. Moreover, since the centroid shift is qualitatively similar to our 
p(p k ) J p( pk ^ anc l the axial ratio is qualitatively related to our P2/P0, it would seem unlikely that 
we would reach entirely opposite conclusions. 

The other plausible variable is to consider how BTb and Mohr et al. computed their statistics 
on the real cluster data. The ROSAT data analyzed by BTb have better spatial resolution and 
sensitivity than the Einstein data analyzed by Mohr et al.. This implies that the Mohr et al. data 
should be biased in the direction of less "substructure" with respect to BTb, which is the opposite 
of what is found. Another important difference between the two investigations is that the PR m 
are computed within apertures of fixed metric size, whereas Mohr et. al. use a S/N criterion to 
define the aperture size. The fixed metric radius used by the PR m ensures that cluster structure 
on the scale ~ R ap is compared consistently which is not true for the S/N criterion (see BTa); 
e.g., Mohr et al. use aperture sizes of 0.38/igo 1 Mpc for Coma and of OMh^ Mpc for A2256. It 



5 This would not necessarily rule out low-density, flat models having Slo < 0.35. 

6 If the small sample of clusters in the Navarro et al. simulation are in fact biased towards more relaxed 
configurations at the present day, then the agreement discussed above between the PR m computed for the Slo = 1 N- 
body simulations in this paper and the PR m that TB computed for the Navarro et al. simulation further strengthens 
the SCDM-ROSAT discrepancy. 
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is not obvious, however, how this confusion of cluster scales would explain the discrepancy of our 
results with Mohr et al..[] 

6. Conclusions 

Using the power ratios (PR m = log 10 (P m /P )) of Buote & Tsai (1995 - BTa; 1996 - BTb; 
Tsai & Buote 1996 - TB) we have examined the sensitivity of galaxy cluster morphologies to 
and P(k) using large, high-resolution N-body simulations. X-ray images are generated from the 
dark matter by letting the gas density trace the dark matter. We argue that the PR m should 
not be seriously biased by this approximation because a real gas in hydrostatic equilibrium with 
potentials of CDM clusters is rounder, but also has a larger core radius, the effects of which 
partially cancel. We also argue that the approximation should be reasonable during mergers 
because of the agreement shown between the evolution of the dark matter and gas found by 
TB who analyzed the hydrodynamical simulation of Navarro et al. (1995a). Finally, The PR m 
generated from the N-body simulations in this paper agree with results from the Navarro et al. 
hydrodynamical simulation (TB). Similar agreement is seen between the results of the N-body 
simulations of Jing et al. (1995) and the hydrodynamical simulations of Mohr et al. (1995). 

From analysis of several variants of the standard Cold Dark Matter model, we have shown 
that the PR m can distinguish between models with different f^o and P(k). Generally, influences 
the means of the PR m distributions such that larger values of On primarily imply larger average 
PR values. The slope of the power spectrum and a§ primarily influence the variances of the PR m \ 
smaller n and a$ generally imply smaller PR m variances. 

For examining On, our analysis indicates that PR3 is the best PR m since its mean is quite 
sensitive to f2n but very insensitive to P(k). It is advantageous also to consider PR2 when a 
cosmological constant is introduced since its means differ for the SCDM and LCDM models by 
~ 3a whereas PR3 only distinguishes the models at the ~ 2a level. The dependence of the mean 
of PR2 on P(k) is not overly serious for studying differences in fin because the differences in means 
due to P(k) are always accompanied by larger differences in the variances; i.e. different means 
but consistent variances mostly reflect differences in Qq f° r PRi- Typically, the best apertures for 
segregating models are (0.5, 0.75, 1.0)hg Mpc. 

We did not find it advantageous to compare the distributions along and perpendicular to 
the "evolutionary tracks" in the (PR2, PR4) and (PR2, PR3) planes (see BTb and TB). The 
distributions along the tracks performed essentially as a weighted sum of the constituent PR m - 
The distributions perpendicular to the tracks were in almost all cases consistent for the models. 
Hence, although the evolutionary tracks are useful for categorizing the dynamical states of clusters, 



7 This issue could be addressed by computing PR m on the Einstein sample of Mohr et al., however such a task is 
beyond the scope of the present paper. 
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they do not allow more interesting constraints on and P{k) to be obtained over the individual 
PR m - The consistency of the distributions perpendicular to the evolutionary tracks seems to be a 
generic feature of the CDM models. 

We compared the PR m of the CDM models to the ROSAT cluster sample of Buote & Tsai 
(1996). We find that the means of the Qq = 0.35 OCDM and ROSAT clusters are consistent, but 
the means of PR3 for the LCDM and ROSAT clusters are formally inconsistent at the ~ 3a level. 
We assert that this discrepancy should be considered marginal due to various issues associated 
with the simulation - observation comparison. 

However, the means of PR3 for the SCDM and BCDM models (with f^o = 1) exceed the 
ROSAT means by almost a full decade with a high level of significance (~ 4cr). Though the formal 
significance level of this pjj^.j / X-ray comparison should be considered only an approximation, we 
argue that taking into account the hydrodynamics and cooling will not reconcile a discrepancy this 
large. We conclude that the Clo = 1 CDM models cannot produce the observed PR m of the ROSAT 
clusters, and that the discrepancy in PR m means is due to £Iq being too large. This agrees with 
our conclusions obtained in TB for the small sample of clusters drawn from the hydrodynamic 
simulation of Navarro et al. (1995a). These conclusions are also consistent with other indicators 
of a low value of f^o such as the dynamical analyses of clusters (e.g., Carlberg et al. 1995), the 
large baryon fractions in clusters (e.g., White et al. 1993), and the heating of galactic disks (Toth 
& Ostriker 1992). 

Our conclusions are inconsistent with those of Mohr et al. (1995) who instead concluded 
that their Einstein cluster sample favored Qq = 1, CDM over equivalents of our low-density 
models, OCDM and LCDM. We argue that this type of discrepancy is unlikely due to numerical 
differences between our simulations. We discuss possible differences due to how BTb and Mohr et. 
al. computed their statistics on the real cluster data. 

Large hydrodynamical simulations are necessary to render the comparison to the ROSAT 
data more robust. In addition, the effects of combining data at different redshifts needs to be 
explored since cluster formation rates should behave differently as a function of z in different 
models (e.g., Richstone et. al. 1992). It may also prove useful to apply PR m to mass maps of 
clusters obtained from weak lensing (Kaiser & Squires 1993)0, though for cosmological purposes it 
is not clear whether p ma ss will be as responsive as (? gas to different f^o an d P{k). 

We gratefully acknowledge J. Tsai for his role in facilitating this collaboration. We thank 
E. Bertschinger for suggesting to us the scheme to normalize the scale-free models using M*, T. 
Beers for providing his ROSTAT programs to compute Al and TI, and A. Edge for providing the 



s See Wilson, Cole, & Frenk (1996), who have recently studied weak-lensing maps obtained from N-body 
simulations, and concluded that a "global quadrupole statistic" (~ *J P2/P0) can distinguish between low-density 
and critical density models. 
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Table 1: Cosmological Models 



Name 


^0 


Ao 


n 


0"8 


h 


Zi 


SCDM 


1 





1 


1.00 


0.5 


20 


OCDM 


0.35 





1 


0.79 


0.7 


25 


LCDM 


0.35 


0.65 


1 


0.83 


0.7 


39 


BCDM 


1 





1 


0.51 


0.5 


20 


SFOO 


1 














SF10 


1 





-1.0 








SF15 


1 





-1.5 








SF20 


1 





-2.0 









Note. — Zi is the redshift where the simulations started. The scale-free models (SF) are normalized to have the 
same value of M* as SCDM (see §0). 



Table 2: Average Power Ratios 



PRl 



PR* 





0.5 Mpc 


0.75 Mpc 


1.0 Mpc 


0.5 Mpc 


0.75 Mpc 


1.0 Mpc 




avg 


a 


avg 


a 


avg 


a 


avg 


a 


avg 


a 


avg 


a 


SCDM 


-5.14 


0.86 


-5.16 


0.83 


-5.38 


0.76 


-6.72 


0.73 


-6.82 


1.01 


-7.06 


0.98 


OCDM 


-5.55 


0.61 


-5.82 


0.78 


-5.93 


1.03 


-7.40 


0.81 


-7.56 


1.13 


-7.59 


1.28 


LCDM 


-5.45 


0.83 


-5.69 


0.85 


-5.87 


0.83 


-7.08 


0.83 


-7.32 


1.00 


-7.38 


1.02 


BCDM 


-5.01 


0.57 


-5.24 


0.60 


-5.41 


0.63 


-6.83 


0.65 


-6.98 


0.67 


-7.07 


0.68 


SFOO 


-5.34 


0.82 


-5.64 


1.00 


-5.75 


0.97 


-7.12 


1.02 


-7.28 


1.14 


-7.35 


1.10 


SF10 


-5.02 


0.57 


-5.22 


0.71 


-5.47 


0.85 


-6.89 


0.87 


-7.09 


1.02 


-7.36 


1.18 


SF15 


-5.20 


0.85 


-5.22 


0.94 


-5.40 


0.91 


-6.92 


0.99 


-6.93 


1.02 


-7.05 


0.96 


SF20 


-5.07 


0.60 


-5.24 


0.62 


-5.45 


0.57 


-7.01 


0.71 


-7.13 


0.75 


-7.20 


0.79 



Note. - 



Aperture sizes assume h=0.8. 
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Table 3: Significance Tests for Power Ratios 



Models 


0.5 Mpc 0.75 Mpc 1.0 Mpc 

Pt PF Pks Pt PF PKS Pt PF PKS 

(%) (%) (%) | (%) (%) (%) | (%) (%) (%) 


PR 2 


SCDM vs. OCDM 
SCDM vs. LCDM 
OCDM vs. LCDM 
SCDM vs. BCDM 
SFOO vs. SF20 


1.63 4.09 0.94 
10.47 87.90 21.79 
54.64 5.79 51.42 
43.06 1.51 98.09 
10.76 6.35 7.31 


0.06 71.71 0.07 
0.74 89.05 1.97 
47.52 61.72 21.79 
65.19 4.65 21.79 
4.19 0.38 0.94 


0.89 7.18 0.94 
0.83 58.78 3.90 
78.26 20.53 34.56 
83.52 24.87 34.56 
10.44 0.13 3.90 


PR 3 


SCDM vs. OCDM 
SCDM vs. LCDM 
OCDM vs. LCDM 
SCDM vs. BCDM 
SFOO vs. SF20 


0.02 52.99 0.18 
4.44 43.91 21.79 
9.10 88.39 21.79 
49.58 47.26 70.85 
57.65 2.89 21.79 


0.31 46.71 0.94 
3.26 97.10 7.31 
30.99 44.52 51.42 
40.62 1.39 34.56 
49.89 1.28 21.79 


4.43 10.64 7.31 
15.88 81.58 51.42 
43.51 16.61 51.42 
97.46 2.48 70.85 
50.92 4.31 12.97 



Note. — Aperture sizes assume h=0.8. 



Table 4: PR Statistics for ROSAT Clusters 





0.5 Mpc 


1.0 Mpc 




avg a 


avg a 


PR 2 


-5.70 0.44 


-6.00 0.50 


PRs 


-7.62 0.77 


-7.61 0.77 



Note. - 



Aperture sizes assume h=0.8. 
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Table 5: Significance Tests for ROSAT Clusters 



0.5 Mpc 1.0 Mpc 



Models pt(%) p F {%) Pks(%) 


Pt(%) PF {%) Pks(%) 


PR 2 


SCDM 0.60E-01 0.12E-01 0.68E-03 
BCDM 0.10E-04 0.12E+02 0.56E-04 
OCDM 0.23E+02 0.52E+01 0.16E+02 
LCDM O.llE+02 0.21E-01 0.13E+01 


0.14E-01 0.24E+01 0.33E-01 
0.12E-01 0.20E+02 0.12E-01 
0.68E+02 0.23E-01 0.43E+02 
0.41E+02 0.69E+00 0.23E+02 


PR 3 


SCDM 0.15E-03 0.71E+02 0.69E-01 
BCDM 0.64E-03 0.28E+02 0.32E-02 
OCDM 0.22E+02 0.80E+02 0.42E+02 
LCDM 0.44E+00 0.69E+02 0.97E+00 


0.17E+01 0.20E+02 0.10E+02 
0.35E+00 0.46E+02 0.12E+02 
0.93E+02 0.83E+00 0.29E+02 
0.33E+02 0.14E+02 0.32E+02 



Note. - 



Aperture sizes assume h=0.8. 
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Fig. 1.- 




a laa zaa laa saa laa zaa laa saa 

Contour plots of the "X-ray Images" for 16 of the 39 clusters analyzed in the SCDM model 
obtained from projecting j gas oc p 2 DM - Each image is 4 x 4 h~ 2 Mpc 2 and the axes units are 20h 
Mpc pixels. 
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Fig. 4.— 

Joint PR m distributions in the (0.5, 1.0) /igQ 1 
BCDM models. 
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Fig. 5 — 

Histograms for the PRs in the 0.5h^ Mpc aperture. SCDM is given by the solid line, OCDM by 
the dotted line, and LCDM by the dashed line. 




Fig. 6.— 

As Figure ||, but for the 1.0/igQ 1 Mpc aperture. 
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Fig. 7.- 

The standard deviation as a function of the average value of the PRs in the 0.75/igQ 1 Mpc aperture 
for the models in §||. The error bars represent la errors estimated from 1000 bootstrap 
resamplings. 
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Fig. 8.- 

Joint PR m distributions in the (0.5, 1.0) /igQ 1 Mpc apertures for the scale-free models: 
(n = 0) denoted by dots and SF20 (—2) denoted by crosses. 
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Fig. 9.— 

Histograms for the PR m in the 0.5/igg 1 Mpc aperture. Spectral index n = (SFOO) is given by the 
solid line and spectral index n = —2 (SF20) by the dotted line. 
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Fig. 10.— 

As Figure ||, but for the 1-0/igQ 1 Mpc aperture. 
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Fig. 11 — 

The standard deviation as a function of the average value of the PRs in the 0.75/igQ 1 Mpc aperture 
for the models in §|j. The error bars represent la errors estimated from 1000 bootstrap 
resamplings. 



b 0.8 - 




0.6 - 



-5.B -5.6 -5.4 

Average PR 2 



1.2 



1 - 



0.8 



0.6 



I 1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 I 



SF00 



SF10 SF15c 



BCDM 



_L 



_L 



-7.6 -7.4 -7.2 -7 

Average PR 3 



-6.8 



-6.6 



-35- 



Fig. 12 — 

Joint PR m distributions in the (0.5, l.O)/^ 1 Mpc apertures for the ROSAT (crosses) and SCDM 
(dots) clusters. 
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Fig. 13.— 

Histograms for the PRs in the O.bh^ Mpc aperture. ROSAT is given by the solid line, SCDM by 
the dotted line, and OCDM by the dashed line. 
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Fig. 14 — 

As Figure but for the 1.0/igg 1 Mpc aperture. 
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Fig. 15 — 

The standard deviation as a function of the average value of the PRs in the 0.5/igQ 1 Mpc aperture 



for the ROSAT clusters and the models discussed in §^1. The error bars represent la errors 
estimated from 1000 bootstrap resamplings. 
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